Ni-Based Molecular Sieves Nanomaterials for Dry Methane Reforming: Role of Porous Structure and Active Sites Distribution on Hydrogen Production

Global warming, driven by greenhouse gases like CH4 and CO2, necessitates efficient catalytic conversion to syngas. Herein, Ni containing different molecular sieve nanomaterials are investigated for dry reforming of methane (DRM). The reduced catalysts are characterized by surface area porosity, X-ray diffraction, Raman infrared spectroscopy, CO2 temperature-programmed desorption techniques, and transmission electron microscopy. The active sites over each molecular sieve remain stable under oxidizing gas CO2 during DRM. The reduced 5Ni/CBV10A catalyst, characterized by the lowest silica–alumina ratio, smallest surface area and pore volume, and narrow 8-ring connecting channels, generated the maximum number of active sites on its outer surface. In contrast, the reduced-5Ni/CBV3024E catalyst, with the highest silica–alumina ratio, more than double the surface area and pore volume, 12-ring sinusoidal porous channels, and smallest Ni crystallite, produced the highest H2 output (44%) after 300 min of operation at 700 °C, with a CH4:CO2 = 1:1, P = 1 atom, gas hour space velocity (GHSV) = 42 L gcat−1 h−1. This performance was achieved despite having 25% fewer initial active sites, suggesting that a larger fraction of these sites is stabilized within the pore channels, leading to sustained catalytic activity. Using central composite design and response surface methodology, we successfully optimized the process conditions for the 5Ni/CBV3024E catalyst. The optimized conditions yielded a desirable H2 to CO ratio of 1.00, with a H2 yield of 91.92% and a CO yield of 89.16%, indicating high efficiency in gas production. The experimental results closely aligned with the predicted values, demonstrating the effectiveness of the optimization approach.


Introduction
The climate emergency requires significant efforts for effective mitigation.Greenhouse gases (GHGs), such as carbon dioxide and methane, contribute to 90% of total emissions, exacerbating global warming [1].Global temperature increases lead to biodiversity loss, extreme weather events and glacier melting [2,3].Industrial activities have significantly contributed to GHGs, aggravating climate change.As a result, the research on converting CH 4 and CO 2 into H 2 and CO (syngas) has garnered significant global interest.This process, known as the dry reforming of methane (DRM), produces syngas, which has significant applications in both synthetic processes and energy production [4].
Cost-effective Ni and Co as active sites are preferred over costly noble metal active sites; whereas, in the mean of CH 4 interaction energy, Ni is preferred over Co [5,6].Catalytic materials at the nanoscale often have unique catalytic properties to the bulk [7], like strong adsorption energy for CH 4 and CO 2 , easy CO 2 dissociation, and excellent carbon elimination [7].Maria et al. observed that upon moving from bulk to nano, the infrared vibration frequency shifted to a higher magnitude (500 cm −1 to 570 cm −1 ), indicating enhanced metal-support interaction at the nano level [8].The detailed reaction scheme of DRM and other parallel reactions is shown in Figure 1.The DRM reaction can be summed up by two steps, namely dissociation of CH 4 over edge-step sites of Ni and oxidation of dissociated methane (CH 4 -x; x = 1-4) into syngas by CO 2 [9,10] (Figure 1A,B).The DRM reaction is highly endothermic and takes place within 550-850 • C. Within this temperature range, various parallel competitive reactions also run over the catalyst surface which can alter the final H 2 yield and CO yield ratio over the catalyst.Carbon monoxide disproportionation reaction (Boudouard reaction; 2CO → CO 2 +C) is a side reaction and is responsible for carbon deposit.However, by adjusting a temperature of ~700 • C, this reaction becomes thermodynamically infeasible [11] (Figure 1C).The interaction of H 2 (one of the DRM products) with CO 2 cannot be neglected (Figure 1D).Especially, a reverse water gas shift reaction (H 2 +CO 2 → H 2 O +CO) is thermodynamically feasible, which gives water and CO.The consumption of H 2 and the formation of CO in RWGS reactions affect the H 2 /CO ratio seriously.However, the gasification of carbon by water (C + H 2 O → CO + H 2 ) further improves the H 2 and CO yield during the DRM reaction (Figure 1E).
Cost-effective Ni and Co as active sites are preferred over costly noble metal active sites; whereas, in the mean of CH4 interaction energy, Ni is preferred over Co [5,6].Catalytic materials at the nanoscale often have unique catalytic properties to the bulk [7], like strong adsorption energy for CH4 and CO2, easy CO2 dissociation, and excellent carbon elimination [7].Maria et al. observed that upon moving from bulk to nano, the infrared vibration frequency shifted to a higher magnitude (500 cm −1 to 570 cm −1 ), indicating enhanced metal-support interaction at the nano level [8].The detailed reaction scheme of DRM and other parallel reactions is shown in Figure 1.The DRM reaction can be summed up by two steps, namely dissociation of CH4 over edge-step sites of Ni and oxidation of dissociated methane (CH4-x; x = 1-4) into syngas by CO2 [9,10] (Figure 1A,B).The DRM reaction is highly endothermic and takes place within 550-850 °C.Within this temperature range, various parallel competitive reactions also run over the catalyst surface which can alter the final H2 yield and CO yield ratio over the catalyst.Carbon monoxide disproportionation reaction (Boudouard reaction; 2CO → CO2 +C) is a side reaction and is responsible for carbon deposit.However, by adjusting a temperature of ~700 °C, this reaction becomes thermodynamically infeasible [11] (Figure 1C).The interaction of H2 (one of the DRM products) with CO2 cannot be neglected (Figure 1D).Especially, a reverse water gas shift reaction (H2 +CO2 → H2O +CO) is thermodynamically feasible, which gives water and CO.The consumption of H2 and the formation of CO in RWGS reactions affect the H2/CO ratio seriously.However, the gasification of carbon by water (C + H2O → CO + H2) further improves the H2 and CO yield during the DRM reaction (Figure 1E).Alumina is typical support for Ni active sites, but it has limitations for DRM due to its high acidity and the tendency of Ni to migrate into the alumina.Zeolites, which are aluminosilicates, offer a better alternative because they have a rigid microporous framework that supports Ni active sites effectively [12,13].The unique arrangement of pores and windows in zeolites facilitates the diffusion of small molecules, including reactants and products.Ni/USY catalysts show higher reducibility, while Ni/ZSM-5 has a stronger metal-support interaction compared to Ni/Al2O3 [14-17], leading to higher activity and better stability.Synthetic aluminosilicate molecular sieves, which mimic natural zeolites, have tunable pore sizes based on their silica-alumina ratio, allowing for optimization through specific synthetic strategies [18,19].
In addition to using powerful catalysts, achieving optimal catalytic activity in DRM involves fine-tuning various operating parameters such as temperature, the CH4/CO2 Alumina is typical support for Ni active sites, but it has limitations for DRM due to its high acidity and the tendency of Ni to migrate into the alumina.Zeolites, which are aluminosilicates, offer a better alternative because they have a rigid microporous framework that supports Ni active sites effectively [12,13].The unique arrangement of pores and windows in zeolites facilitates the diffusion of small molecules, including reactants and products.Ni/USY catalysts show higher reducibility, while Ni/ZSM-5 has a stronger metal-support interaction compared to Ni/Al 2 O 3 [14-17], leading to higher activity and better stability.Synthetic aluminosilicate molecular sieves, which mimic natural zeolites, have tunable pore sizes based on their silica-alumina ratio, allowing for optimization through specific synthetic strategies [18,19].
In addition to using powerful catalysts, achieving optimal catalytic activity in DRM involves fine-tuning various operating parameters such as temperature, the CH 4 /CO 2 ratio, and space velocity.However, finding the best combination of these parameters through experimentation is labor-intensive, time-consuming, and energy-demanding.Thankfully, statistical tools like response surface methodology (RSM) can predict the optimal catalytic performance using minimal experimental data.RSM continuously exhibits a high level of predictive power, with modeled and experimental results closely matching.RSM includes several designs, such as the Box-Behnken design (BBD), central composite design (CCD), and Daehler design (DD) [20].For instance, RSM predicted a 30.98%H 2 yield and 69.36% CO yield for an yttrium oxide-cobalt/mesoporous alumina catalyst at 900 • C reaction temperature, 30,690 mLgcat −1 h −1 WHSV, and 1.027 CH 4 /CO 2 [21].Experimental results under the same conditions showed yields of 32.12% H 2 and 67.88% CO.Similarly, for calcium iron oxide-supported nickel catalyst, RSM predicted 77.82% H 2 yield and 75.76%CO yield at 832.45 • C reaction temperature, 0.96 CH 4 /CO 2 ratio, and 35,000 mL g −1 h −1 gas hour space velocity [22], closely matching experimental results.
These close correlations between predictions and experiments were also observed for other catalysts, including Ta-promoted Ni/ZSM-5, alumina-magnesia supported Ni-W, and dendritic fibrous SBA-15-supported Ni catalysts, using the central composite design under RSM [23][24][25].Today, RSM is not limited to optimizing the reaction parameters but also to fine-tuning the catalyst's composition, such as the proportion of active sites, promoters, and supports, to achieve the best catalytic activity for DRM reactions [26].
The work aims to produce a new, practical, and dependable Ni-based catalyst for the DRM process supported by synthetic molecular sieves with a distinct pore structure and a variable silica-alumina ratio.Unlike conventional techniques, this novel approach seeks to improve catalytic performance and stability.Several sophisticated instrumentation techniques are used to confirm the catalytic activity results, such as X-ray diffraction, surface area porosity, Raman infrared spectroscopy, and temperature-programmed reduction/desorption. Furthermore, we utilized an advanced statistical method, response surface methodology (RSM), to systematically investigate the simultaneous effects of various operating parameters on the optimized catalyst's performance.Finally, we thoroughly evaluate the catalyst's long-term stability to ensure sustained and efficient production of H 2 and CO, underscoring the practical applicability and robustness of our cutting-edge catalyst design.

Catalyst Preparation
The synthetic molecular sieves (CBV3024E, CP810B, CBV 10A, CBV 20A) were made using a traditional wet impregnation (WI) technique; each had a 5 wt.% metal loading of nickel nitrate hexahydrate (Ni (NO 3 ) 2 •6H 2 O, Alfa Aesar; Thermo Fisher Scientific, Lancashire, UK) catalysts.The calcined molecular sieves were gradually treated with a water-based solution of metal precursors.The slurry was made by adding support to the solution while stirring at 80 • C for 2 h, followed by drying at 120 • C overnight, and then calcining for three hours at 600 • C. The catalysts are identified as follows: 5Ni/CBV 3024E, 5Ni/CP810B, 5Ni/CBV 10A, and 5Ni/CBV 20A.

Catalyst Performance Evaluation
An apparatus made by PID Eng & Tech, with dimensions of 9.1 mm in diameter and 30 cm in length, is used to conduct the DRM reaction experiment with 0.1 g of catalyst.
Encircling the reactor is a furnace that generates heat.Located axially in the catalyst bed's center, a K-type thermocouple measures the bed's temperature.The catalyst must be reduced for one hour at 600 • C in a stream of flowing hydrogen gas (at a flow rate of 30 mL/min) before it can be used.Using a combination of CH 4 , CO 2 , and N 2 at flow rates of 30 mL/min, 30 mL/min, and 10 mL/min, the DRM reaction was conducted over a reduced catalyst at 700 • C and atmosphere pressure of 1.An online gas chromatography system, which includes a thermal conductivity detector (TCD), Porapak Q columns, and a molecular sieve 5A, is used to examine the gas streams entering and leaving the system.The following is the formula for determining the H 2 yield and CO yield:

Systematic Testing for Improvement
In this study, central composite design (CCD) (Stat-Ease Design Expert software, version 13) was employed with three independent variables that included reaction temperature (A), CH 4 /CO 2 ratio (B), and space velocity (SV) (C) to predict the optimum process conditions for the dry reforming of CH 4 over the optimum-obtained 5Ni/CBV3024E during the 300 min on stream.Fifteen trials were conducted, three times each, as part of the experimental design and optimization.

Catalyst Characterization
The Rigaku (Miniflex) diffractometer, operating at 40 mA and 40 kV, was used to perform the powder X-ray diffraction (XRD) investigation of the newly prepared catalyst using Cu Kα radiation (λ = 0.15406 nm).A laser Raman spectrometer from JASCO in Japan was used for the Raman analysis.It had a 532 nm beam excitation and 1.6 mW laser intensity.The sample was exposed for 10 s at 3 accumulations.Quantification of graphitization and carbon deposition type on catalysts was completed with the use of a laser Raman spectrometer (NMR-4500-JASCO, Tokyo, Japan) manufactured by JASCO.A 532 nm wavelength excitation laser was utilized.The IR Prestige-21 SHMADZU was used to record the Fourier transform infrared (FTIR) spectra.Surface characteristics and functional groupings of the samples were investigated using the measures.To conduct the H 2 -TPR, a 0.07 g sample was heated to 1000 • C at a rate of 10 • C/min while being cycled through a gas mixture of H 2 /Ar (v/v, 10/90) at a rate of 40 mL/min.A thermal conductivity detector was used to monitor the signal of H 2 consumption.Carbon dioxide (CO 2 ) temperature-programmed desorption was used to examine the decreased catalyst's basicity using the same detector.At room temperature, a mixture of 10% CO 2 /He gas was allowed to interact with 70 mg of catalyst.Various strengths of CO 2 adsorption over the catalyst were based on the strengths of the basic sites.Depending on the CO 2 binding strength to the surface, desorbed CO 2 temperatures were different when the temperature was raised to 1000 • C. Using TCD, the amount of CO 2 that was adsorbed was analyzed.This cycle of H 2 TPR-CO 2 TPD-H 2 TPR persisted constantly.

Determinations of Porosity Parameters
The specific surface area (SSA), pore volume, and average pore diameter were determined using nitrogen adsorption-desorption isotherms at liquid N2 temperature (77 K) on a Micromeritics TriStar II plus instrument.The BET equation P/ ] was employed to calculate SSA, where P = the equilibrium pressure of the gas; P0 = is the saturation pressure of the gas; V = is the volume of gas adsorbed; Vm = is the volume of gas corresponding to a complete monolayer cover-age of the adsorbent; C = is a constant related to the energy of adsorption.Pore volume was calculated by subtracting the monolayer volume from the total adsorbed gas volume at a relative pressure of approximately 0.99.The average pore diameter was derived from the total pore volume and specific surface area using the relation Dpore = 4 Vpore/S.The pore size distribution was determined using the BJH method from the desorption branch of the isotherm.Before analysis, samples were degassed at 300 • C under vacuum for 24 h to remove any adsorbed contaminants.

Catalyst Name
The Crystallite Size of Ni (nm) 5Ni/CBV 3024E 13.9 5Ni/CBV810B 13.9 5Ni/CBV10A 15.9 5Ni/CBV20A 22.2 The surface area and porosity result of reduced-5Ni/CBV3024E, reduced-5Ni/CP810B, reduced-5Ni/CBV10A, and reduced-5Ni/CBV20A catalysts are shown in Figure 3 and Table 2.The adsorption isotherm of all catalysts belongs to type IV, having an H3 hysteresis loop, which is characterized by infinite adsorption at a high relative pressure (p/p o ), indicating the presence of slit-like pores [30].The surface area of the Ni-incorporated molecular sieve is significantly less than the pristine molecular sieve, indicating the deposition of Ni crystallite into the pores of the sieve.The forced closure along the desorption branch is observed at different relative pressures for each catalyst [31].This is known as the tensile strength effect.The closure point at a certain relative pressure is due to the desorption of adsorbate from pores at which stress in the adsorbate reaches tensile strength [32].Reduced-5Ni/CBV10A catalyst belongs to a mordenite porous framework having a minimum SiO2/Al2O3 ratio (13), minimum surface area (142 m 2 /g), and pore volume (0.06 cm 3 /g).Such inferior porous architecture in 5Ni/CBV10A results in an early closure point of the desorption curve (at 0.04 p/p o ).5Ni/CBV20A catalyst again belongs to The surface area and porosity result of reduced-5Ni/CBV3024E, reduced-5Ni/CP810B, reduced-5Ni/CBV10A, and reduced-5Ni/CBV20A catalysts are shown in Figure 3 and Table 2.The adsorption isotherm of all catalysts belongs to type IV, having an H3 hysteresis loop, which is characterized by infinite adsorption at a high relative pressure (p/p°), indicating the presence of slit-like pores [30].The surface area of the Ni-incorporated molecular sieve is significantly less than the pristine molecular sieve, indicating the deposition of Ni crystallite into the pores of the sieve.The forced closure along the desorption branch is observed at different relative pressures for each catalyst [31].This is known as the tensile strength effect.The closure point at a certain relative pressure is due to the desorption of adsorbate from pores at which stress in the adsorbate reaches tensile strength [32].Reduced-5Ni/CBV10A catalyst belongs to a mordenite porous framework having a minimum SiO 2 /Al 2 O 3 ratio (13), minimum surface area (142 m 2 /g), and pore volume (0.06 cm 3 /g).Such inferior porous architecture in 5Ni/CBV10A results in an early closure point of the desorption curve (at 0.04 p/p°).5Ni/CBV20A catalyst again belongs to the mordenite framework, having a higher SiO 2 /Al 2 O 3 ratio (20).The surface parameter of reduced-5Ni/CBV20A is improved two times more (390 m 2 /g surface area and 0.13 cm 3 /g pore volume) than reduced-5Ni/CBV10A.Reduced-5Ni/CP810B catalyst has a ZSM-11-type porous framework, having a higher SiO 2 /Al 2 O 3 ratio ( 25), but the most enriched surface parameters like surface area (496 m 2 /g) and pore volume (0.62 cm 3 /g).Reduced-5Ni/CBV3024E catalyst comprises the highest ratio of silica/alumina (30), and its pore architect belongs to the ZSM-8 framework.It attains two times the surface area (304 m 2 /g) and pore volume (0.13 cm 3 /g) of the reduced-5Ni/CBV10A catalyst.The closure point for reduced-5Ni/CBV20A and reduced-5Ni/CBV3024 catalyst falls between 0.41 and 0.42 p/p°.The pore size distribution plot (dV/dlog(w) vs. W; where "V" and "w" are pore volume and pore width) shows the presence of multimodal pores from the microporous-to-macroporous range (Figure S2).However, the average pore diameter of supports and catalysts vary from 7.52 to 16.58 nm and 5.55 to 15.70 nm in the mesoporous range, respectively.pore volume) than reduced-5Ni/CBV10A.Reduced-5Ni/CP810B catalyst has a ZSM-11type porous framework, having a higher SiO2/Al2O3 ratio (25), but the most enriched surface parameters like surface area (496 m 2 /g) and pore volume (0.62 cm 3 /g).Reduced-5Ni/CBV3024E catalyst comprises the highest ratio of silica/alumina (30), and its pore architect belongs to the ZSM-8 framework.It attains two times the surface area (304 m 2 /g) and pore volume (0.13 cm 3 /g) of the reduced-5Ni/CBV10A catalyst.The closure point for reduced-5Ni/CBV20A and reduced-5Ni/CBV3024 catalyst falls between 0.41 and 0.42 p/p o .The pore size distribution plot (dV/dlog(w) vs. W; where "V" and "w" are pore volume and pore width) shows the presence of multimodal pores from the microporous-tomacroporous range (Figure S2).However, the average pore diameter of supports and catalysts vary from 7.52 to 16.58 nm and 5.55 to 15.70 nm in the mesoporous range, respectively.The Raman spectra of reduced-5Ni/CBV 3024E, reduced-5Ni/CBV810B, reduced-5Ni/CBV 10A, and reduced-5Ni/CBV 20A catalysts are shown in Figure 4. Interestingly, upon reduction, the intensity of the Raman spectra of reduced-5Ni/CBV3024E catalyst is increased prominently, whereas the rest of the catalysts' selected vibration modes of Raman spectra are slightly intensified upon reduction (Figure S3).That means the polarizability pattern over the catalyst is modified greatly upon reduction under hydrogen and it is specified as the aluminosilicates being the architect of molecular sieve.The Raman vibration bands corresponding to the bending and starching vibration modes of O-Si-O and Si-O have been reliably identified at approximately 511 cm −1 [33] and 854 cm −1 [34].Cubic alumina is itself Raman, inactive [35], but it perturbs the stretching vibrations of silicates under the alumino-silicate framework, which results in a series of Raman bands for vibration of Si-O at 283 cm −1 , 403 cm −1 , 923 cm −1 , 1051 cm −1 , and 1166 cm −1 and vibration bands for Al-O at 679 cm −1 [36,37].The Raman vibrational frequency at 139 cm −1 is attributed to the translation of Si/Si and Al/Al lattice modes [35,38].The additional Raman bands at 446 cm −1 and 553 cm −1 are attributed to ring vibration in the zeolite framework [39].The reduced-5Ni/CBV 3024E catalyst has the highest Si/Al ratio so it shows the most intense Raman bands.In the same way, reduced-5Ni/CBV10A catalyst has a minimum Si/Al ratio and so it has minimum Raman band intensity for Si-O and Al-O vibration.The other two catalysts, reduced-5Ni/CP810B and reduced-5Ni/CBV20A, have comparable Raman bands (Figure 4B).increased prominently, whereas the rest of the catalysts' selected vibration modes of Raman spectra are slightly intensified upon reduction (Figure S3).That means the polarizability pattern over the catalyst is modified greatly upon reduction under hydrogen and it is specified as the aluminosilicates being the architect of molecular sieve.The Raman vibration bands corresponding to the bending and starching vibration modes of O-Si-O and Si-O have been reliably identified at approximately 511 cm −1 [33] and 854 cm −1 [34].Cubic alumina is itself Raman, inactive [35], but it perturbs the stretching vibrations of silicates under the alumino-silicate framework, which results in a series of Raman bands for vibration of Si-O at 283 cm −1 , 403 cm −1 , 923 cm −1 , 1051 cm −1 , and 1166 cm −1 and vibration bands for Al-O at 679 cm −1 [36,37].The Raman vibrational frequency at 139 cm −1 is attributed to the translation of Si/Si and Al/Al lattice modes [35,38].The additional Raman bands at 446 cm −1 and 553 cm −1 are attributed to ring vibration in the zeolite framework [39].The reduced-5Ni/CBV 3024E catalyst has the highest Si/Al ratio so it shows the most intense Raman bands.In the same way, reduced-5Ni/CBV10A catalyst has a minimum Si/Al ratio and so it has minimum Raman band intensity for Si-O and Al-O vibration.The other two catalysts, reduced-5Ni/CP810B and reduced-5Ni/CBV20A, have comparable Raman bands (Figure 4B).The infrared spectra of reduced-5Ni/CBV3024E, reduced-5Ni/CP810B, reduced-5Ni/CBV10A, and reduced-5Ni/CBV20A catalysts are shown in Figure 5.The molecular sieve CBV3024E is of the ZSM type, which may be identified by the presence of an infrared band with a wavelength of approximately 548 cm −1 (Band A) and 448 cm −1 (Band B) [40].The intensity ratio of the peak at 548 cm −1 (IA) and 448 cm −1 (IB) for the reduced-5Ni/CBV3024E catalyst is found to be 0.8.The IA/IB > 0.6 is characterized as ZSM-8 zeolite [28] (Figure 5A).In the reduced-5Ni/CP810B catalyst, the peak at 548 cm −1 is split into two peaks at 528 cm −1 and 572 cm −1 .Here, the vibration peak intensity at 528 cm −1 is selected as Band A. The IA/IB ratio for the reduced-5Ni/CP810B catalyst is found to be 0.28.IA/IB < 6 is the characteristic ratio for ZSM-11 among the ZSM family.Reduced-5Ni/CBV10A and reduced-5Ni/CBV20A exhibit typical vibration patterns of mordenite at 460 cm −1 , 565 cm −1 , 800 cm −1 , and 1100 cm −1 [41].The infrared spectra of reduced-5Ni/CBV3024E, reduced-5Ni/CP810B, reduced-5Ni/ CBV10A, and reduced-5Ni/CBV20A catalysts are shown in Figure 5.The molecular sieve CBV3024E is of the ZSM type, which may be identified by the presence of an infrared band with a wavelength of approximately 548 cm −1 (Band A) and 448 cm −1 (Band B) [40].The intensity ratio of the peak at 548 cm −1 (I A ) and 448 cm −1 (I B ) for the reduced-5Ni/CBV3024E catalyst is found to be 0.8.The I A /I B > 0.6 is characterized as ZSM-8 zeolite [28] (Figure 5A).In the reduced-5Ni/CP810B catalyst, the peak at 548 cm −1 is split into two peaks at 528 cm −1 and 572 cm −1 .Here, the vibration peak intensity at 528 cm −1 is selected as Band A. The I A /I B ratio for the reduced-5Ni/CP810B catalyst is found to be 0.28.I A /I B < 6 is the characteristic ratio for ZSM-11 among the ZSM family.Reduced-5Ni/CBV10A and reduced-5Ni/CBV20A exhibit typical vibration patterns of mordenite at 460 cm −1 , 565 cm −1 , 800 cm −1 , and 1100 cm −1 [41].The cyclic H2TPR-CO2TPD-H2TPR results of different Ni-containing molecular sieves are shown in Figure 6 and Table 3.The reduction peak below 300 °C is found to be negative due to hydrogen spillover into pores of Ni-containing molecular sieves [42,43].The H2- The cyclic H 2 TPR-CO 2 TPD-H 2 TPR results of different Ni-containing molecular sieves are shown in Figure 6 and Table 3.The reduction peak below 300 • C is found to be negative due to hydrogen spillover into pores of Ni-containing molecular sieves [42,43].The H 2 -TPR profile above 300 • C is due to the consumption of H 2 by reducible metal oxide over the catalysts.In H 2 -TPR, the 5Ni/CBV10A catalyst absorbs the maximum amount of H 2 (20 cm 3 /g), indicating the presence of the maximum amount of reducible NiO in the catalyst.Interestingly, H 2 is consumed mainly at about 350 • C and a little at high temperatures.This peak is attributed to the reduction in Ni 2+ , which is secluded within the charge compensation sites of zeolite [44,45].The DRM reaction is carried out over a reduced catalyst, and the reduced catalyst may contain basicity for possible the interaction of acidic CO 2 gas during the DRM reaction.The CO 2 -TPD profile of the reduced catalyst showed that most of the CO 2 desorption occurs (2.2 cm 3 /g) at about 100 • C, which is typically attributed to weak basic sites or basicity contributed by surface hydroxyl [46].At about 275-400 • C, there is again desorption of some CO 2 (0.71 cm 3 /g), which indicates that basic sites of moderate strength (basicity contributed by surface oxide) are also present over the catalyst surface [47][48][49].Interestingly, after CO2 desorption of the reduced catalyst, if H2-TPR is again carried out, there is no reduction peak.It indicates that CO2 is an oxidizing gas, and it can oxidize the carbon deposit during the DRM, but it will not oxidize the catalytic active sites (metallic Ni).That means active sites remain stable under the oxidizing gas environment during DRM.
The H2 consumption over 5Ni/CBV20A and 5Ni/CP810B are comparable (15-15.7  Interestingly, after CO 2 desorption of the reduced catalyst, if H 2 -TPR is again carried out, there is no reduction peak.It indicates that CO 2 is an oxidizing gas, and it can oxidize the carbon deposit during the DRM, but it will not oxidize the catalytic active sites (metallic Ni).That means active sites remain stable under the oxidizing gas environment during DRM.
The H 2 consumption over 5Ni/CBV20A and 5Ni/CP810B are comparable (15-15.7 cm 3 /g) and it is about 25% less than 5Ni/CBV10A.So, the amount of reducible NiO species is proportionally less over 5Ni/CBV20A and 5Ni/CP810B catalysts.Among 5Ni/CBV20A and 5Ni/CP810B, a relatively higher amount of H 2 is consumed at a high-temperature region (600 • C) over 5Ni/CP810B.The reduction peak of NiO in zeolite at about 500 • C was reported for such NiO particles, which were located inside the zeolite channel [44,45,50].So, here, the reduction peak of about 600 • C can be attributed to the reduction in strong interactions inside the pore channels of the molecular sieve.Overall, it can be said that the amount of reducible species is similar in 5Ni/CP810 and 5Ni/CBV20A, but the earlier one has a greater amount of NiO inside the pore channel of the molecular sieve.The CO 2 -TPD profile of these reduced catalysts shows depletion of low-temperature desorption peaks (concerning 5Ni/CBV10 catalyst), indicating the absence of weak basic sites due to surface hydroxyl.However, the basic sites of moderate strength are present over both catalysts considerably.After CO 2 -TPD treatment of reduced catalyst, if H 2 -TPR is again run, there is no reducible peak.It indicates that the active sites remain stable in the presence of oxidizing gas like CO 2 .The reducibility pattern of 5Ni/CBV3024E is similar to 5Ni/CP810B but it is shifted more towards a relatively lower reduction temperature.That means the 5Ni/CBV3024E catalyst has a higher edge of reducibility than the 5Ni/CP810B catalyst.The basicity pattern of the reduced-5Ni/CBV3024E catalyst is similar to the reduced-5Ni/CP810B catalyst.Here also, the active sites remain stable against oxidizing gas CO 2 (as observed in the rest of the catalysts).
The grainy and featureless textures of 5Ni/CBV3024E, 5Ni/CBV20A, and 5Ni/CP810B are evident in SEM images (Figure S4).The TEM image and particle size distribution of reduced and spent 5Ni/CBV3024E catalysts are shown in Figure 7. Analysis of the reduced catalysts using X-ray diffraction (XRD) confirms the conversion of NiO to metallic Ni.TEM images show larger, darker features that could be indicative of larger Ni particles on the catalyst surface.Smaller, fainter contrasts might correspond to smaller particles within the zeolite pores.However, due to TEM contrast limitations, definitive location assignment is difficult.The average particle size of Ni is found at 6.47 nm over a reduced-5Ni/CBV3024E catalyst.In the spent-5Ni/CBV3024E, the average particle size grows to 8.44 nm.Interestingly, no carbon nanostructures like nanofiber and nano-tubes are observed over the catalyst that was used.the catalyst surface.Smaller, fainter contrasts might correspond to smaller particles within the zeolite pores.However, due to TEM contrast limitations, definitive location assignment is difficult.The average particle size of Ni is found at 6.47 nm over a reduced-5Ni/CBV3024E catalyst.In the spent-5Ni/CBV3024E, the average particle size grows to 8.44 nm.Interestingly, no carbon nanostructures like nanofiber and nano-tubes are observed over the catalyst that was used.

Discussion
Ni containing different molecular sieves are investigated for DRM reforming of methane.XRD, Raman, and IR spectra of reduced catalyst samples confirm that CBV10A and CBV20A belong to mordenite-based support, whereas CP810B and CBV3024E belong to ZSM-11 and ZSM-8-based supports, respectively.The silica-alumina ratios of CBV10A, CBV20A, CP810B, and CBV3024E are 13, 20, 25, and 30, respectively.ZSM-5, ZSM-11, and modernities are categorized under pore channel types.The porous channel of mordenite is made up of 12-ring channels and 8-ring channels where 12-ring channels are connected through 8-ring channels arranged in parallel and perpendicular fashion [51].ZSM-11 is weaved by intersectional straight 12-ring channels, whereas ZSM-8 has intersectional straight and sinusoidal 12-ring channels [52].

Discussion
Ni containing different molecular sieves are investigated for DRM reforming of methane.XRD, Raman, and IR spectra of reduced catalyst samples confirm that CBV10A and CBV20A belong to mordenite-based support, whereas CP810B and CBV3024E belong to ZSM-11 and ZSM-8-based supports, respectively.The silica-alumina ratios of CBV10A, CBV20A, CP810B, and CBV3024E are 13, 20, 25, and 30, respectively.ZSM-5, ZSM-11, and modernities are categorized under pore channel types.The porous channel of mordenite is made up of 12-ring channels and 8-ring channels where 12-ring channels are connected through 8-ring channels arranged in parallel and perpendicular fashion [51].ZSM-11 is weaved by intersectional straight 12-ring channels, whereas ZSM-8 has intersectional straight and sinusoidal 12-ring channels [52].
CBV20A and CP810B contain mordenite and ZSM-11-type pores channels.Interestingly, upon impregnating 5 wt.%Ni over this molecular sieve, NiO is dispersed over the outer surface of the molecular sieve (at charge compensation sites) and inside the pore channel of the molecular sieves (confirmed by H2 TPR and TEM).Approximately the same amount of reducible NiO is reduced into metallic Ni over these molecular sieves throughout the reduction in catalysts.However, in 5Ni/CP810B, there are more active sites incorporated into the pore channel than the 5Ni/CBV20A catalyst.Clearly, this difference is CBV20A and CP810B contain mordenite and ZSM-11-type pores channels.Interestingly, upon impregnating 5 wt.%Ni over this molecular sieve, NiO is dispersed over the outer surface of the molecular sieve (at charge compensation sites) and inside the pore channel of the molecular sieves (confirmed by H 2 TPR and TEM).Approximately the same amount of reducible NiO is reduced into metallic Ni over these molecular sieves throughout the reduction in catalysts.However, in 5Ni/CP810B, there are more active sites incorporated into the pore channel than the 5Ni/CBV20A catalyst.Clearly, this difference is due to more dispersion of smaller Ni crystallites (13.9 nm) into the 12-ring channels of the 5Ni/CP810B catalyst (ZSM-11 architect).In the case of reduced-5Ni/CBV20A (SiO 2 /Al 2 O 3 = 20), porous channels are the mordenite type where 12-ring channels are connected through 8-ring channels.The crystallite size of Ni is also the largest at 22.2 nm.So, passing large Ni crystallites into "8-ring channels" or "12-ring channels through 8-ring channels" is relatively difficult.Overall, the reduced-5Ni/CBV20A catalyst has a smaller number of active sites in the channel, whereas the reduced-5Ni/CP810B catalyst has a larger population of active sites in the channel.
Additionally, the reduced-5Ni/CP810B catalyst has the highest surface area and pore volume among all catalysts.The active sites inside the channel are quite stable compared to the active sites on the outer surface of the molecular sieve.Reduced-5Ni/CBV20A catalyst, which has less amount of stable active site into mordenite's porous architects, shows 27% H 2 yield constantly and 36-34% CO yield during 300 min on stream.However, reduced-5Ni/CP810B, which has a higher population of stable active sites in the ZSM-11's porous architects, performs better.Initially, 41% H 2 yield and 50% CO yield were observed over reduced-5Ni/CP810B, which slowed down to 37% H 2 yield and 46% CO yield during 300 min on stream, respectively.The presence of higher CO yield than H 2 yield over both catalysts confirms the existence of hydrogen-consuming reactions like reverse water gas shift reactions as major side reactions (Figure 8C) [50].
5Ni/CBV10A catalyst also has a mordenite-based porous framework but a lower silica-alumina ratio (SiO 2 /Al 2 O 3 = 13).Increasing the proportion of alumina with silica was reported to reduce the size of NiO aggregates, resulting in enhanced reducibility and high dispersion [53].XRD results also show a smaller Ni crystallite (15.9 nm) over reduced-5Ni/CBV10A than 22 nm Ni crystallite over reduced-5Ni/CBV20A.It indicates a higher dispersion of Ni crystallites over a reduced-5Ni/CBV10A catalyst.5Ni/CBV10A catalyst has 33% more reducible NiO-species than 5Ni/CBV20A and 5Ni/CP810B catalysts.The minimum surface area and pore volume, early closure point of desorption, and narrower 8-ring channels of the mordenite framework make localization of Ni difficult into the pore channels.The H 2 -TPR result justifies the presence of most of the reducible species at the outer surface of the molecular sieve which is stabilized at the charge compensation sites of the molecular sieve.Despite the deficit of active sites in the pore channels, the reduced-5Ni/CBV10A catalyst acquires the highest amount of total active sites and achieves 42% H 2 yield and 53% CO yield during the 300 min time on stream.The porous architect of reduced-5Ni/CBV3024E and reduced-5Ni/CP810B are like ZSM-8 and ZSM-11 types.The former has straight intersectional sinusoidal 12-ring channels, but the latter has simply straight intersectional 12-ring channels.The Ni's crystallinity, reducibility profile, and basicity profile of reduced-5Ni/CBV3024E is very similar to reduced-5Ni/CP810B.But the former has a higher edge of reducibility than the latter.Overall, the catalytic performance of reduced-5Ni/CBV3024E is best and it acquires 44% H 2 yield and 52% CO yield after up to 300 min on stream.Now, the reaction scheme for DRM over the efficient Ni-supported molecular sieve, one having minimum SiO 2 /Al 2 O 3 ratio ( 13) and another having maximum SiO 2 /Al 2 O 3 ratio (30), can be illustrated (Figure 9). Figure S5 shows the TGA profiles of the used samples after 5 h of reaction.The figure displays that the 5Ni/CBV10A sample exhibits the greatest carbon deposition, whereas the carbon deposit is the least over 5Ni/CBV20A as well as 5NiCBV3024E catalysts.Reduced-5Ni/CBV10A (SiO2/Al2O3 = 13) has mordenite-type pores architecture, and it acquires the least surface area and pore volume.Over this catalyst, localization of metallic Ni into 12-ring channels connected through narrow 8-ring channels is quite difficult due to space constraints and minimum pore volume.However, the maximum population of active sites has been created at the outer surface of the molecular sieve (Figure 9A), which catalyzes the decomposition of the C-H bond (of CH4).The sequential oxidation of decomposed CH4 by CO2 conveys 42% hydrogen yield after up to 300 min.Reduced-5Ni/CBV3024E (SiO2/Al2O3 = 13) has a ZSM-8-type porous architecture (made up of sinusoidal 12-ring channels) and the smallest Ni crystallites resulting in localization of Ni into the pores mostly (Figure 9B).Despite the low concentration of active sites, the active sites clipped under pore channels of reduced-5Ni/CBV3024E catalyst are quite stable.The active sites (the outer surface as well as inside the pore channels) catalyze C-H decomposition, which further undergoes oxidation under CO2.In the TEM image of this catalyst, carbon nanostructures are also not observed.The catalyst achieves 44% H2 yield in up to 300 min.This investigation validates the remarkable performance of the 5Ni/CBV3024E catalyst, building on earlier studies (mentioned in Table S1).Reduced-5Ni/CBV10A (SiO 2 /Al 2 O 3 = 13) has mordenite-type pores architecture, and it acquires the least surface area and pore volume.Over this catalyst, localization of metallic Ni into 12-ring channels connected through narrow 8-ring channels is quite difficult due to space constraints and minimum pore volume.However, the maximum population of active sites has been created at the outer surface of the molecular sieve (Figure 9A), which catalyzes the decomposition of the C-H bond (of CH 4 ).The sequential oxidation of decomposed CH 4 by CO 2 conveys 42% hydrogen yield after up to 300 min.Reduced-5Ni/CBV3024E (SiO 2 /Al 2 O 3 = 13) has a ZSM-8-type porous architecture (made up of sinusoidal 12-ring channels) and the smallest Ni crystallites resulting in localization of Ni into the pores mostly (Figure 9B).Despite the low concentration of active sites, the active sites clipped under pore channels of reduced-5Ni/CBV3024E catalyst are quite stable.The active sites (the outer surface as well as inside the pore channels) catalyze C-H decomposition, which further undergoes oxidation under CO 2 .In the TEM image of this catalyst, carbon nanostructures are also not observed.The catalyst achieves 44% H 2 yield in up to 300 min.This investigation validates the remarkable performance of the 5Ni/CBV3024E catalyst, building on earlier studies (mentioned in Table S1).

Central Composite Design (CCD)
CCD is a frequently utilized approach in response surface methodology (RSM), aimed at extracting the maximum amount of information about a process from a limited number of experiments.This technique is effective in modeling quadratic effects and determining optimal conditions for process enhancement based on the desired properties of the design and the number of factors considered.This study's experimental variables include temperature, the CH 4 :CO 2 ratio, and space velocity, each with specified upper and lower limits.The mean value of x i is (X oi ) and is given in Equation ( 4), which describes the center of the interval of the lower (x i min ) and upper limits (x i max ), and the deviation of either limit from the mean is described in Equation ( 5): In experimental design, the coordinates (X o1 , X o2 , X o3 , . .., X on ) define the center point of the experiment.This particular set of input variables is generally incorporated to detect curvature in the response surface.The initial values are converted into coded, dimensionless variables Xi within the range of −1 to 1 using the transformation specified in Equation ( 6).
where x i , represent the original value of the i-th input variable, and ∆x i signifies the range of values for that variable.
In experimental design, the outcome variable is often expressed as a polynomial equation of several factors.The Taylor Series expansion is used to derive this polynomial equation, approximating the response surface with a mix of linear, quadratic, and occasionally cubic terms.In our experiment, a modified quadratic polynomial is found, and the full quadratic model with three factors in the general form is given by Equation ( 7): where, X 1 , X 2 , and X 3 are the experimental factors in actual or coded values, β 0 is a constant term, β i , i = 1, 2, 3 are the linear terms, β ii are the quadratic terms, β ij , j = 1, 2, 3 are the interaction terms, and ε is the error term, which represents the random variation in Ŷ that is not accounted for by the process parameters [39].Table 4 lists the actual and coded values of the studied factors.

Process Modeling and Analysis of Variance
ANOVA, or analysis of variance, is a statistical technique that partitions the total variability in data into distinct components attributed to various sources.This method is critical in identifying significant factors and their interactions, enhancing a good theoretical model.Process modeling and variance analysis utilize the central composite design (CCD) approach, which optimizes processes, assesses factor interactions, and determines the optimal power transformation for response data to normalize or stabilize variance [40].Analysis of variance for the various components is shown in Table 5. High F-values and low p-values indicate significant model terms at a confidence level of approximately 95%.Moreover, high R 2 values indicate that the model effectively represents experimental data, as illustrated in Figure 10.detect curvature in the response surface.The initial values are converted into coded, dimensionless variables Xi within the range of −1 to 1 using the transformation specified in Equation ( 6).
where x , represent the original value of the i-th input variable, and ∆x signifies the range of values for that variable.
In experimental design, the outcome variable is often expressed as a polynomial equation of several factors.The Taylor Series expansion is used to derive this polynomial equation, approximating the response surface with a mix of linear, quadratic, and occasionally cubic terms.In our experiment, a modified quadratic polynomial is found, and the full quadratic model with three factors in the general form is given by Equation ( 7): where, X , X , and X are the experimental factors in actual or coded values, β is a constant term, β , i = 1,2,3 are the linear terms, β are the quadratic terms, β , j = 1,2,3 are the interaction terms, and ε is the error term, which represents the random variation in  that is not accounted for by the process parameters [39].Table 4 lists the actual and coded values of the studied factors.ANOVA, or analysis of variance, is a statistical technique that partitions the total variability in data into distinct components attributed to various sources.This method is critical in identifying significant factors and their interactions, enhancing a good theoretical model.Process modeling and variance analysis utilize the central composite design (CCD) approach, which optimizes processes, assesses factor interactions, and determines the optimal power transformation for response data to normalize or stabilize variance [40].Analysis of variance for the various components is shown in Table 5. High F-values and low pvalues indicate significant model terms at a confidence level of approximately 95%.Moreover, high R 2 values indicate that the model effectively represents experimental data, as illustrated in Figure 10.

Final Equations and Accuracy Models
According to the experimental data, ANOVA was conducted at a significance level α = 0.05, and after identifying the significant factor effects while excluding insignificant ones, the following models were developed using Design-Expert software version 13.
H 2 yield = −554.42115+ 1.86263A − 172.05000B − 0.004409 C + 0.301000 AB + 0.0000051 AC − 0.001332 A 2 − 59.60000 B 2 (8) CO yield = 117.66724− 0.031819 A − 3.81000 B − 0.003599 C + 0.0000044 AC H 2 /CO = 1.20391 − 0.001375 A − 0.280110 B + 0.0000018 C + 0.002200 AB + −0.0000054 BC − 0.514286 B 2 (10) where A is temperature, B is the feed ratio of reactants, and C is the SV, which is the hourly space velocity.Within the models ( 5)-( 7), the constant coefficients signify the anticipated value of the response variable when all other factors are at zero.The coefficient of one factor indicates its main effect, illustrating how the response variable changes with a one-unit increase in each factor while holding all other factors constant.coefficient associated with the interaction of two factors represents the combined effects of these interactions.Positive coefficients indicate a positive or increased influence on the response variable, while negative coefficients indicate the opposite effect.
The coefficient of determination (R 2 ) quantifies the proportion of variability in the response variable that is accounted for by the independent variables in the model, serving as a measure of the model's goodness-of-fit.A high R 2 value indicates a well-fitting model.The R 2 values for the expected H 2 /CO, H 2 yield, and CO yield models are 0.9923, 0.9795, and 0.8846, respectively.These values indicate that the models explain approximately 99.23%, 97.95%, and 88.46% of the total variations in the expected responses.
Table 6 presents the experimental and predicted values of all response variables based on the proposed models.As depicted in Table 6, the proximity of the estimated responses derived from the proposed models to the actual values highlights the models' capacity to accurately reflect the observed data.This is evidenced by the low absolute error rates of the three models, indicating that, on average, predicted values deviate from actual values by 0.78%, 1.04%, and 1.59%, respectively, as shown in Table 6.A lower Mean Absolute Percentage Error (MAPE) value signifies greater accuracy in the models.In addition, MAPE, Absolute Error Percentage (APE), and Mean Absolute Error (MAE) which are shown in Equations ( 8)- (10) are used to evaluate the model's accuracy.Table 6 confirms a strong correlation between CCD models and experimental results, with R 2 near 1.This consistency is evident in Figure 10; through plotting predicted against actual values that a crucial for model assessment, we found the close alignment to the X = Y line, indicating a good fit.
where APE = Absolute Percentage Error; E i : Actual observed value; P i : Predicted value from the model.
MAE = Mean Absolute Error; n = Number of data points.
MAPE = Mean Absolute Percentage Error.

Simulation of Design-Expert Program
Due to the endothermic nature of the DRM reaction, higher temperatures significantly boosted the yields of both H 2 and CO.This can be explained by the increased kinetic energy within the system at higher temperatures.This rise in energy effectively accelerates the movement of reactants (mass transfer), leading to a more efficient reaction and, ultimately, greater production of the desired products.
It was observed that a decrease in product yield alongside an increase in the CH 4 /CO 2 ratio from 0.5 to 1.0.This suggests that excess CH 4 is decomposing and depositing carbon, leading to lower overall product output.Additionally, the study found that yields of both H 2 and CO decreased with a higher reactant ratio and increased space velocity.Interestingly, the ratio of H 2 to CO increased with a higher reactant ratio but decreased with a higher space velocity.

One Factor Effect (2D) Plot
The effect of each process parameter on the reaction responses is shown in Figures 11-13.Figure 11 indicates that increasing temperature, increasing the ratio, and decreasing the SV value will increase H 2 /CO. Figure 12 indicates that increasing temperature, decreasing the ratio, and decreasing the SV value will increase H 2 yield.Figure 13 indicates that increasing temperature, decreasing the ratio, and decreasing the SV value will increase CO yield.Due to the endothermic nature of the DRM reaction, higher temperatures significantly boosted the yields of both H2 and CO.This can be explained by the increased kinetic energy within the system at higher temperatures.This rise in energy effectively accelerates the movement of reactants (mass transfer), leading to a more efficient reaction and, ultimately, greater production of the desired products.
It was observed that a decrease in product yield alongside an increase in the CH4/CO2 ratio from 0.5 to 1.0.This suggests that excess CH4 is decomposing and depositing carbon, leading to lower overall product output.Additionally, the study found that yields of both H2 and CO decreased with a higher reactant ratio and increased space velocity.Interestingly, the ratio of H2 to CO increased with a higher reactant ratio but decreased with a higher space velocity.

One Factor Effect (2D) Plot
The effect of each process parameter on the reaction responses is shown in Figures 11-13.Figure 11 indicates that increasing temperature, increasing the ratio, and decreasing the SV value will increase H2/CO.Figure 12 indicates that increasing temperature, decreasing the ratio, and decreasing the SV value will increase H2 yield. Figure 13 indicates that increasing temperature, decreasing the ratio, and decreasing the SV value will increase CO yield.Two Factors Effect (3D Plot) Employing the regression coefficients for various components led to the development of a series of response surface methodology (RSM) equations.With the assistance of these equations and the Design-Expert software, response surface plots were generated to predict the conversion or formation of components in the reaction system across two varying process variables while holding the third variable constant.These predictions are illustrated in the 3D models shown in Figures 14-19.Two Factors Effect (3D Plot) Employing the regression coefficients for various components led to the development of a series of response surface methodology (RSM) equations.With the assistance of these equations and the Design-Expert software, response surface plots were generated to predict the conversion or formation of components in the reaction system across two varying process variables while holding the third variable constant.These predictions are illustrated in the 3D models shown in Figures 14-19.Two Factors Effect (3D Plot) Employing the regression coefficients for various components led to the development of a series of response surface methodology (RSM) equations.With the assistance of these equations and the Design-Expert software, response surface plots were generated to predict the conversion or formation of components in the reaction system across two varying process variables while holding the third variable constant.These predictions are illustrated in the 3D models shown in Figures 14-19.Figures 14 and 15 show the three-dimensional response surface plot, which sents the effects of the factors (temperature, SV, and ratio CH : CO on the varia H2/CO.Figure 14 shows the surface plots which represent the relationship betw response variable (H2/CO) and the two factors (temperature and the ratio CH : CO = 24000.It is shown that by increasing the temperature and increasing the ratio CH the H2/CO increases.Figure 15 shows the surface plots that represent the function tionship between a designated response variable (H2/CO) and the two-factor va (Temperature and SV) with CH : CO fixed at 0.89.The response surface shows th increasing the temperature and decreasing the SV, the H2/CO increases.It was ob to increase from 0.67 at 750 °C to 1 at 850 °C.
Figures 16 and 17 show the three-dimensional response surface plots that rep the effects of the factors (temperature, SV, and ratio CH : CO ) on the variation in th of H2. Figure 16 shows the surface plots which represent the relationship between sponse variable (H2 yield) and the two factors, temperature and the ratio CH : CO 24000.It is shown that by increasing the temperature and decreasing the ratio, H increases. Figure 17 shows the surface plots that represent the functional relations tween a designated response variable (H2 yield) and the two factor variables, temp and SV, with CH : CO fixed at 0.8081.The response surface shows that by increas temperature and decreasing the SV, the yield of H2 increases.It was observed to i from 69.31% at 750 °C to 95.85% at 850 °C.
Figures 18 and 19 show the three-dimensional response surface plots that rep the effects of the factors (temperature, SV, and ratio CH : CO ) on the variation yields.Figure 18 shows the surface plots which represent the relationship betw Figure 14 shows the surface plots which represent the relationship between the response variable (H 2 /CO) and the two factors (temperature and the ratio CH 4 :CO 2 ) at SV = 24,000.It is shown that by increasing the temperature and increasing the ratio CH 4 :CO 2 , the H 2 /CO increases.Figure 15 shows the surface plots that represent the functional relationship between a designated response variable (H 2 /CO) and the two-factor variables (Temperature and SV) with CH 4 :CO 2 fixed at 0.89.The response surface shows that with increasing the temperature and decreasing the SV, the H 2 /CO increases.It was observed to increase from 0.67 at 750 • C to 1 at 850 • C.
Figures 16 and 17 show the three-dimensional response surface plots that represents the effects of the factors (temperature, SV, and ratio CH 4 : CO 2 ) on the variation in the yield of H 2 .Figure 16 shows the surface plots which represent the relationship between the response variable (H 2 yield) and the two factors, temperature and the ratio CH 4 :CO 2 , at SV = 24,000.It is shown that by increasing the temperature and decreasing the ratio, H 2 yield increases. Figure 17 shows the surface plots that represent the functional relationship between a designated response variable (H 2 yield) and the two factor variables, temperature and SV, with CH 4 :CO 2 fixed at 0.8081.The response surface shows that by increasing the temperature and decreasing the SV, the yield of H 2 increases.It was observed to increase from 69.31% at 750 • C to 95.85% at 850 • C.
Figures 18 and 19 show the three-dimensional response surface plots that represents the effects of the factors (temperature, SV, and ratio CH 4 :CO 2 ) on the variation in CO yields.Figure 18 shows the surface plots which represent the relationship between the response variable (CO yield) and the two factors, temperature and the ratio CH 4 :CO 2 , at SV = 24,000.
It is shown that with increasing temperature and decreasing ratio, the CO yield increases. Figure 19 shows the surface plots that represent the relationship between the response variable CO yield and the two factors variables, temperature and SV, with CH 4 :CO 2 fixed at 0.6209.The plot shows that with increasing temperature and decreasing SV, the CO yield increases.It was observed to increase from 74.09% at 750 • C to 92.15% at 850 • C.

Conclusions
Ni containing different molecular sieves, named CBV10A, CBV20A, CP810B, and CBV3024, are investigated for DRM reaction.Reduced-5Ni/CBV10A catalyst has the lowest silica-alumina ratio and a mordenite-type porous framework, exhibiting a unique behavior.Due to its minimal pore volume, early desorption closure point, and narrow pore channels (8-ring), the Ni active sites are primarily localized on the outer surface, achieving a 42% hydrogen yield over 300 min of continuous operation.
In contrast, the other catalysts, namely, reduced-5Ni/CBV20A, reduced-5Ni/CP810B, and reduced-5Ni/CBV3024E, displayed 25% fewer total active sites, which were distributed both on the outer surface and within the pore channels.Despite reduced-5Ni/CBV20A having double the surface area and pore volume compared to reduced-5Ni/CBV10A, its larger Ni crystallites and narrower 8-ring channels restricted the localization of active sites within the pores.This resulted in a lower hydrogen yield of 27%, as the total number of active sites was inferior.
The reduced-5Ni/CP810B catalyst, featuring broader 12-ring channels of the ZSM-11 type and the smallest Ni crystallites, allowed for more stable active sites within the pores, achieving a better hydrogen yield of 41%.Similarly, the reduced-5Ni/CBV3024E catalyst, with its sinusoidal 12-ring channels and highest silica-alumina ratio, exhibited similar reducibility and crystallinity profiles to reduced-5Ni/CP810B.However, the active sites on reduced-5Ni/CBV3024E formed at relatively lower temperatures, leading to the highest hydrogen yield of 44% over 300 min.
To optimize the process parameters, we used a simulation of the Design-Expert program based on the CCD to develop a quadratic model for RSM and to optimize the process parameters.The ANOVA table obtained suggested that the model was significant.A set of 15 experiments was carried out to study optimization at different reaction conditions.The findings revealed that the reaction temperature had the greatest impact on the process, followed by the space velocity, while the feed ratio had little impact.It was discovered that the ideal parameters for maximizing the response variables were a reaction temperature of 850 • C, a space velocity of 22,000.009ccg −1 h −1 , and a feed ratio of 0.932.Under these conditions, the projected H 2 and CO yields were 91.5% and 89.37%, respectively, and the H 2 /CO ratio was 1.00.The experimental results obtained under the predicted optimal operating conditions matched the predicted values quite well, and hence, the validation experiment confirmed the model's accuracy with negligible differences between predicted and observed results.

Figure 1 .
Figure 1.The detailed reaction scheme of DRM and other parallel reactions over supported Ni-based catalyst.

Figure 1 .
Figure 1.The detailed reaction scheme of DRM and other parallel reactions over supported Ni-based catalyst.

Figure 10 .
Figure 10.The experimental and the predicted data for (a) H2/CO, (b) H2 yield, and (c) CO yield.Figure 10.The experimental and the predicted data for (a) H 2 /CO, (b) H 2 yield, and (c) CO yield.

Figure 10 .
Figure 10.The experimental and the predicted data for (a) H2/CO, (b) H2 yield, and (c) CO yield.Figure 10.The experimental and the predicted data for (a) H 2 /CO, (b) H 2 yield, and (c) CO yield.

Figure 11 .
Figure 11.The relationship between the reaction parameters and the H2/CO ratios.Figure 11.The relationship between the reaction parameters and the H 2 /CO ratios.

Figure 11 .
Figure 11.The relationship between the reaction parameters and the H2/CO ratios.Figure 11.The relationship between the reaction parameters and the H 2 /CO ratios.

Figure 11 .
Figure 11.The relationship between the reaction parameters and H2/CO ratios.

Figure 12 .
Figure 12.The relationship between the reaction parameters and H2 yield percentage.Figure 12.The relationship between the reaction parameters and H 2 yield percentage.

Figure 12 . 27 Figure 13 .
Figure 12.The relationship between the reaction parameters and H2 yield percentage.Figure 12.The relationship between the reaction parameters and H 2 yield percentage. Nanomaterials 2024, 14, x FOR PEER REVIEW 19 of 27

Figure 14 .
Figure 14.The relationship between the temperature and H2/CO ratio at SV = 24,000.

Figure 13 .
Figure 13.The relationship between the reaction parameters and CO yield percentage.

Nanomaterials 2024 , 27 Figure 13 .
Figure 13.The relationship between the reaction parameters and CO yield percentage.

Figure 14 .
Figure 14.The relationship between the temperature and H2/CO ratio at SV = 24,000.Figure 14.The relationship between the temperature and H 2 /CO ratio at SV = 24,000.

Figure 14 .
Figure 14.The relationship between the temperature and H2/CO ratio at SV = 24,000.Figure 14.The relationship between the temperature and H 2 /CO ratio at SV = 24,000.

Figure 14 .
Figure14.The relationship between the temperature and H2/CO ratio at SV = 24,000.

Figure 15 .
Figure 15.The relationship between the temperature, SV, and H2/CO at ratio = 0.89.

Figure 18 .
Figure18.The relationship between the temperature, ratio, and CO yield % at SV = 24,000.

Figure 19 .
Figure 19.The relationship between the temperature, SV, and CO yield % at Ratio = 0.6209.Figures 14 and 15 show the three-dimensional response surface plot, which represents the effects of the factors (temperature, SV, and ratio CH 4 :CO 2 on the variation in H 2 /CO.Figure14shows the surface plots which represent the relationship between the response variable (H 2 /CO) and the two factors (temperature and the ratio CH 4 :CO 2 ) at SV = 24,000.It is shown that by increasing the temperature and increasing the ratio CH 4 :CO 2 , the H 2 /CO increases.Figure15shows the surface plots that represent the functional relationship between a designated response variable (H 2 /CO) and the two-factor variables (Temperature and SV) with CH 4 :CO 2 fixed at 0.89.The response surface shows that with increasing the temperature and decreasing the SV, the H 2 /CO increases.It was observed to increase from 0.67 at 750 • C to 1 at 850 • C.Figures16 and 17show the three-dimensional response surface plots that represents the effects of the factors (temperature, SV, and ratio CH 4 : CO 2 ) on the variation in the yield of H 2 .Figure16shows the surface plots which represent the relationship between the response variable (H 2 yield) and the two factors, temperature and the ratio CH 4 :CO 2 , at SV = 24,000.It is shown that by increasing the temperature and decreasing the ratio, H 2 yield increases.Figure17shows the surface plots that represent the functional relationship between a designated response variable (H 2 yield) and the two factor variables, temperature and SV, with CH 4 :CO 2 fixed at 0.8081.The response surface shows that by increasing the temperature and decreasing the SV, the yield of H 2 increases.It was observed to increase from 69.31% at 750 • C to 95.85% at 850 • C.Figures18 and 19show the three-dimensional response surface plots that represents the effects of the factors (temperature, SV, and ratio CH 4 :CO 2 ) on the variation in CO yields.Figure18shows the surface plots which represent the relationship between the response variable (CO yield) and the two factors, temperature and the ratio CH 4 :CO 2 , at SV = 24,000.

Figure 20 .
Figure 20.Optimum predicted simultaneous values of H2/CO ratio, H2 yield, and CO yield at the controlled conditions in Table7.

Figure 20 .
Figure 20.Optimum predicted simultaneous values of H 2 /CO ratio, H 2 yield, and CO yield at the controlled conditions in Table7.

Author
Contributions: A.S.A.-F.: Project administration, methodology, formal analysis, review, data curation, writing-original draft, editing.A.A.I.: Methodology, conceptualization, investigation, formal analysis, writing-original draft.A.H.F.: Investigation, resources, supervision, review and editing.A.I.O.: Conceptualization, funding acquisition, resources, formal analysis, data curation, writing-review and editing.Y.M.A.: Formal analysis, data curation, visualization.F.S.A.: Validation, analysis, writing-review and editing.A.E.A.: Methodology, software, data curation, investigation.All authors have read and agreed to the published version of the manuscript.Funding: Deputyship for Research and Innovation, "Ministry of Education" in Saudi Arabia.Data Availability Statement: Data are contained within the article.

Table 1 .
The crystallite sizes of the different zeolite-supported catalysts.

Table 1 .
The crystallite sizes of the different zeolite-supported catalysts.

Table 2 .
The surface area, pore volume, and average pore size of different reduced catalysts and their corresponding supports.

Table 2 .
The surface area, pore volume, and average pore size of different reduced catalysts and their corresponding supports.

Table 3 .
The amount of H 2 consumption and CO 2 desorption against different temperatures under H 2 TPR-CO 2 TPD-H 2 TPR cyclic experiments.

Table 4 .
Actual and coded values for the process parameters.

Table 5 .
Analysis of variance for the various components.

Response 3: CO yield (R 2 = 0.8846)
Where A represents the temperature, B is the CH 4 /O 2 ratio, and C represents the space velocity.

Table 4 .
Actual and coded values for the process parameters.

Table 6 .
Experimental and predicted data results for various components of the reaction system.

Table 7 .
Comparison of theoretical model predictions and experimental findings.

Table 7 .
Comparison of theoretical model predictions and experimental findings.

Table 8 .
Comparison of catalytic activity of other catalysts reported in the previous literature.

Table 8 .
Comparison of catalytic activity of other catalysts reported in the previous literature.